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ABSTRACT 


To meet the increasing energy needs and to optimize the use 
of fossil fuel, new sources of energy were developed. However all 
of these sources are not safe and easy to use. Out of these 

hydrogen is an important source of energy. It can be stored in 
large amount in metal hydrides to transport and use safely. The 
basic property by which metal hydrides store hydrogen and release 
it on heating is of importance. Enthalpy of formation play an 
important role in stability of these hydrides and intermetallic 
compounds. Only those hydrides are stable for which AH is 
negative. Partial replacement of metal hydrides with metal like 
( V / Cr / Mn / Ni ) gives pronounced effect on the stability. 

In present work , enthalpy of formation of various compounds 
and hydrides of Ti-Fe-H system is calculated using Fortran 
Programming and effect of partial replacement of Fe by ( V / Ni / 

Ti / Cr ) on stability has been studied theoretically. 

The absorption and desorption of H by intermetallic is 

St 

associated with release & absorption of heat energy. One can think 
of using therefore, metal hydrides to pump heat from one location 
to other easily. An attempt therefore has been made to study the 
feasibility of indigenously developed Mischmetal based 

intermetallic for heat pumping application. 
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CHAPTER 1 


INTRODUCTION : 

These days, in environment conscious world, there is hue and 
cry to save the environment, control the pollution and optimize 
the use of conventional sources of energy in order to save the 
future generation from suffering of harmful effects. With the 
development, need of energy is increasing day by day. Things which 
were considered earlier as luxury have become necessity today. 
There is need to develop new energy sources in order to 
accommodate the rising energy needs. Keeping with this line^ 
scientist of the world have developed nuclear energy and non- 
conven tiona l energy sources like wind, tidal, solar, wave etc. But 
these sources have not yet become very popular because techniques 
to use these sources are either very expensive or not very 
convenient. But with time, when cost of conventional energy 
increases new sources will become competitive. Further to save the 
earth from pollution and ozone layer depletion new sources will 
become popular. 

Unlike oil and gas, the alternative non-conventional sources 
are site specific. Energy generated from them need to be 
transmitted and distributed after converting into suitable form. 
In general, primary thermal and mechanical energy is converted 
into electricity and transmitted to the user-centers. But for some 
applications such as aviation this method is not suited at all. 
Ideally the secondary fuel should have total use cycle which is 
economic, energy efficient, safe and environmentally acceptable. 
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In addition it should be 

(i) derivable from abundant and renewable resources 

(ii) easily storable and movable over long distances 

(iii) safer in wider range of applications. 

The alternative fuels are either .liquid hydrocarbons, 

» methanol, ethanol and hydrogen ^ 3 . Of these, synthetic 

liquid hydrocarbons are very expensive because of heat cost and 

coal requirement. It’s only advantage lies in it’s offering 
continuity of existing petrofuel distribution system. But this is 
only temporary advantage. Nevertheless it is open to the same 
objection as natural petroleum with regard to environment 
pollution. The production of methane, methanol and ethanol is also 
dependent on the availability of raw material and unfavorable 
energy balance. In this view it is doubtful if any one of these 

can be produced on a large scale to match the usage of the 

petroleum. 

In contrast hydrogen possess several unique advantages over 
the other alternative and satisfies practically all the criterion 
for ideal secondary fuels enumerated earlier. Hydrogen is a 
promising medium for both, energy transmission and storage. It is 
essentially non polluting .Since major byproduct of combustion being 
water (H 0). Further hydrogen has the highest energy density per 
unit weight of any chemical. This can be seen from tablet 1 . 1 ) <4> . 1 1 
has a diversified number of uses ranging from combustion engine to 
fuel cells. One of the greatest limitations in exploitation of 
hydrogen as a fuel however, appears to be difficulty in storing it 
economically and conventionally. Hydrogen can be extracted from 
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TABLE 1.1 

COMPARISION OF THE ENERGY DENSTY IN VARIOUS FUELS 


POUER SOURCE 

ENERGY DENSITY 

CONVERSION 

NET UATT HOUR/KG 


(UATT HOUR/KG) 

EFFICIENCY 


Pb/acid battery 

30 

70 

21.0 

Present 

30 

70 

21.0 

Advanced 

50 

70 

35.0 

Li/MS battery 

150 

70 

105.0 

FeTiH 

i . ? 

516 

30 

154.0 

Mg NiH 

2 * 

1121 

30 

336.0 

MgH 2 C10% Ni) 

2555 

30 

767.0 

Gasoline 

12880 

23 

j 

2962.0 

Gaseous H 

2 

33900 

30 

10170.0 

(100 atm pr . ) 


j 
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fossils fuels and /or water. The most common method of storing 
hydrogen is in gaseous form in cylinders at high pressures. But 
this gives problem of hydrogen embitterment and diffusion 
leading to certain amount of risk. This method is not very safe. 
Other methods, in use to store hydrogen, include absorption by 
metal hydrides. 

Metal hydrides are ■e h'gfrteV afrl compounds of metal and hydrogen. 

Hydrogen is stored safely inside the metal itself in large 
quantities. The amount of hydrogen is greater than either gaseous 
or liquid storage. The major difficulty, however lies in the 
weight and cost of hydriding materials. 

Hydrogen storage density of hydrides often exceeds the 
density of liquid hydrogen. The best example of such hydride is 
water itself. In which partial hydrogen density is 110 Kg/m 3 at 
20°K. The relative partial hydrogen density of some hydrogen 
compounds is shown in table 1.2 tt> . 

Hydride formation reaction is reversible. Also at a given 
temperature, in ideal situation, the absorption (formation of 
hydrides) and desorption (release of hydrogen) take place at 
nearly the same pressure. Practically, however the absorption and 
desorption pressures are not the same at a given temperature. This 
effect is called hysteresis. The temperature for which the 
absorption and desorption take place at atmospheric pressure may 
be very different for different hydrides. The suitability of an 
intermettalic depends on factors such as structure, diffusion rate 
of hydrogen during charge and discharge and enthalpy of hydride 


formation. 



COMPARISON OF HYDROGEN 


COMPOUND 

MOLECULAR 

WEIGHT 

CKfi/Kmole) 

H 0 

2 

18.0 

H SO 

2 4 

98.1 

liqCH^ 

16.0 

liqH z 

2.0 

TiH 

2 

49.9 

LaNi H 

S <5 

438.5 

TiFeH 

i . S>5 

105.7 


5 


LE 1.2 


TS IN VARIOUS HYDROGEN COMPOUNDS 


NO. OF ATOMS 

PARTIAL DENSITY 

PER UNIT VOLUME 

OF HYDROGEN 

%X 0- 2B (m- a ]D 


6.7 

111 

2 . 2 

36 

6.3 

105 

4.2 

71 

9 . 2 

153 

5.3 

88 

6 . 2 

101 
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1 . 1 ) CHARACTER I ST ICS OF METAL HYDRIDE SYSTEM : 

The reaction of hydrogen with metal (M) to make hydride M 
H is mostly exothermic and it can be written as <4> 

X 

m + |h = m h + ah 

The reverse reaction is desorption process and forward 
reaction is absorption process. Reverse reaction is endothermic in 
nature and external heat is to be supplied for desorption 
process. The absorption and desorption process of metal- hydride 
system are studied by means of pressure-composition isotherm in 
fig l.l <4> .It is plot of pressure and composition at various 
temperature . The initial steep slop corresponds to hydrogen going 
into solid solution and this single phase region is a denoted as 
a-phase. The position where the curve begins to change slope onthe 
P-C-T diagram denotes the appearance of a metal hydride of 
/5-phase. Due to varied degree of solubility of hydrogen in 
metals, the resulting metal hydrides sometime become 
non-stoichiometric . In second phase, hydrogen pressure remains 
constant and a plateau results, as more hydrogen is added. The 
plateau continues as long as there are two distinct phases as 
required by Gibb’s phase rule. 

F = C - P + 2 

Where 

F = degree of freedom 

C = number of compositions 

P = number of phases 

The effect of temperature is shown m fig 1.1 -In general, as 

the temperature is increased, the miscibility gap tends to narrow 



HYOROGEN PRESSURE 
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Fig 1.1 


PRESSURE VS HYDROGEN TO METAL RATIO ISOTHERM 
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and gap disappears and o*~phase converts continuously into ft- 
phase. Yet another hydride phase(^) may be formed, for which a 
higher plateau pressure result^ 

Plateau pressure which varies strongly with temperature can 


be described approximately in a limited temperature range by 


<i> 


In p = 

H 

2 


AS AH 
R RT 


Uhere 


AH = change in enthalpy/ mole of hydrogen 
AS = change in entropy/ mole of hydrogen 
R = universal gas constant 

When AS and AH are taken to be temperature independent in the 
limited temperature range considered, a plot of In p vs 1/T (fig 

( 4 ) 

1.2) yields the values for AS and AH. Slope of curve gives A H and 
intercept with 1/T axis gives AS. The negative entropy effect 
predominantly determined by difference between entropy of hydrogen 


in gaseous state and 

in 

hydride (solid state) 

. Since 

the 

later 

energy consideration 

is 

relatively small, -AS 

will 

not 

be 

much 

different from entropy 

contained in the gas 

which 

is 

about 31 


calories/degree mole at one atm. at room temperature. 

The sequence of elementary steps which occur during 
absorption of hydrogen with in a single phase region of solubility 
is 


(1 ) 


1/2 H T ra ?P -°.r_t -JlgLg -g-H Li 5$ e 1/2 H (Surface) 


( 2 ) 

Chemisorption 


■+ H 


(3) 

surface migration 


* H (entering site) 
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Fig 1.2 


VARIATION OF PRESSURE WITH TEMPERATURE 
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(4) 

Chemisorption to 
absorption transitio 

The energy required to break hydrogen molecule into atom is 
equal to 104.20 Kcals/mol H 

2 

Desorption is reverse of the above sequence with the same 
steps as for absorption. Step (3) may not be necessary for some 
sol id /H^ system. The above sequence is limited to reaction 
occurring under isothermal conditions. 

The hydrogen saturated solid and hydride phase (ot and ft) 
co-exist in the plateau region. The sequence of steps can be 
represented as 

1/2 H z — > 1/2 ^(surface adsorption) » H H — \ [H] 

^ ► Interface ^ ► [H] Hydride 

Diffusion (step 5) of hydrogen presumably occurs through the 
hydride to the internal interface where the phase transformation 
is initiated (step 6). Any of these steps may be slow and will 
determine the kinetics of the reaction for a given system. In 
various solubility ranges at different temperatures and pressures 
the rate of determining step may be different as well. 

1.2) HYSTERESIS : 

Many metal hydrides show hysteresis effect (absorption and 
desorption occur at different pressures) at different temperatures 
as shown in P-C-T diagram for absorption-desorption process (fig 
1.3) <4> . The cause of hysteresis is not fully known yet. It ie 


^H] (adjacent to surface) 


Solid 


(5) 

diffusion 


(bulk) 


[H] 
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Fie 1.3 HYSTERESIS PLOT 


OF HYDRIDES 
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believed that increase in lattice parameters (expansion effect) on 
hydr iding play a vital role. An isothermal cycle of 
absorption-desorption is initiated by equi libriat ing gaseous 
hydrogen with the metallic sample at increasing pressures. The 
solid expands due to interstitial occupancy of hydrogen (ot-phase) 
as a function of increasing hydrogen gas pressure. As a result 
micro-regions of the remaining unfilled interstitial sites come in 
a state of increasing compressive strains, as hydrogen proceeds 
until /5-phase starts developing. Thus the oi-phase is subjected to 
plastic deformation at the micro-level due to compressive stresses 
during hydrogen absorption and initiation of /5-phase 
precipitation. The system is henceforth in two phase region now, 
further absorption only increases the amount of /5-phase at cost of 
ot-phase maintaining a constant gas pressure (plateau region). The 
specific volume of /5-phase is considerably lower than ot-phase, 
causing severe macro-strains. The corresponding stresses even 
exceed the fracture strength of bulk sometimes. In general, the 
original material is overall brittle to some extent. The 
brittleness is further enhanced by hydrogenation and the solid 
crumbles into fine powder in the two phase region. The strain is 
high in the /5-phase due to large volume change, that causes 
localised plastic deformation and develops residual compressive 
stresses. After oi-phase is converted into the strained /5-phase 
further hydrogen can be introduced by increasing the gas pressure 
again. The /5-phase matrix is now under increasing compressive 
stress in the micro region, similar to the a-phase where the 
interstitial sites are progressively filled. 
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Hysteresis becomes apparent in desorption as the boundary of 
the two phase (a+/3) region is approached, Dur i ng hydrogen removal 
the /3-phase matrix is increasingly relieved of compressive 
stresses. This causes the desorption isotherm pressure to drop 
below plateau found during the absorption cycle. As one approaches 
the two-phase boundary the /3-phase progr essively get unstrained . The 
onset of ot-phas e precipitation actual ly contributes in the removal 
of strains. The /3-phase progressively disappears in the two-phase 
region. Subsequently the equilibrium pressure decreases as the 
hydrogen is removed from the oi-phase. The cycle is completed when 
hydrogen is removed completely. The degree of hysteresis in a 
metal-hydrogen system is affected such that it is more pronounced 
in dynamic than inconvent ional static tests. Hysteresis reduces the 
efficiency of hydrogen storage system. 

1.3) SURVEY OF HYDRIDES (BINARY AND TERNARY ) <4> : 

The binary compounds of hydrogen with Boron Gallium and 

other elements of the 4 th to 7 th main groups of the periodic table 

are gaseous/volatile. For example CH , NH , HO and HC1. The 

4 a 2 

halogen hydrides when brought into water, form a positive 
univalent ion of hydrogen. The hydrides AlH^, SnH^ and the 
germanium hydrides are typical molecular metal hydride. They have 
a high degree of volatility and a low melting point and are 
thermally unstable. 

The binary compounds of the transition metals, if they 
combine with hydrogen at all, form metallic hydrides. These 
hydrides do not always have composition fixed by stoichiometric 
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ratio. The metals forming such hydrides often first incorporate 
hydrogen purely in solid solution, without change of crystal 
structure, although there is some expansion of the crystal 
lattice. With many metals however, a structural transformation 
sets in when the amount of hydrogen taken up exceeds a certain 
concentration. Since this is associated with a discontinuous 
change in properties it is considered that compound is formed. 
These metallic hydrides are usually very hard and brittle. 

A hydride formed from a binary disordered alloy can be 
considered as a solid solution of two binary hydrides. For such 
cases there is only a small difference from the hydride of single 
metal. A by-hydride formed from an ordered compound, however for 
example LaNi H formed from LaNi , generally has no relation at 

5 o 5 

all to the binary hydrides of the constituents La and Ni . Since 
these hydrides are formed by absorption of hydrogen at room 
temperature, where diffusion of La and Ni atom is negligible, the 
crystal structure of hydride can be derived from that of the 
original compound by relatively small changes in the La and Ni 
position. As a matter of fact the lattice has simply to expand to 
accommodate the hydrogen. The system will show binary behavior 
e.g. 

LaNi + 3H ► LaNi H 

5 2 5 <S 

For compounds that absorb hydrogen at high temperature for 
instance above 300°C, the diffusion of the metal atom may become 
possible leading to ternary behavior e.g. 

3 Zr Cu + 4.75 H = 5 ZrH + ZrCu tt> 

2 2 1. 1> a 
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1.4) CHARACTERISTICS OF SOME IMPORTANT INTERMETTALIC HYDRIDES 

(a) AB Compounds : 

ss 

These hydrides have equilibrium pressure of few atmosphere 5 at 
temperature upto 100°C, low hysteresis and easy activation in the 
initial cycle. However both the metallic constituents of the 
hydrides are very expensive and hence increase the overall cost. 
AB family exhibits hexagonal or orthorhombic structure. LaNi is 

o S5 

most important member of this family. It'S reaction with hydrogen 
is 

LaNi + 3H = LaNi H a> 

5 2 3 <S 

Plateau pressure at room temperature has been found to be 
equal to 2.2 atmosphere. Ni addition decreases the stability of 
rare earth hydrides and improves the capacity of hydrogen 
absorption. To reduce the cost Ni can be replaced by A1 upto some 
extent. Material is then characterized by lower decomposition 
pressure without impairing the kinetics of the hydrogen carrying 
capacity. Replacement of Ni by Pd, Co, Fe, Cr , Ag, Cu , In, Sn and 
Ga bring down the plateau pressure in all cases except for Pd. 
However these substitutions decrease the absorption capacity. 

(b) AB Compounds : 

These intermettalic store hydrogen at the low cost. These 

materials have equilibrium pressure of a few atmospheres at 

o 

temperature upto 100 C. Fe Ti and FeTi are main compounds. Fe Ti 

2 2 

does not absorb appreciable amount of hydrogen. For FeTi to react 
with hydrogen at practical rates reversely, it is necessary to 
activate the material. The activation is accomplished by out 
gassing the alloy at 300°C first and contacting it with hydrogen 
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at a pressure of one atm. later. 

Unlike the AB system the reaction kinetics of AB compound is 

25 

only practically controlled by heat transfer, consequently the 
chemical reaction rates are of much more significance than heat 
transfer capability. 

The Fe-Ti-H system exhibits pronounced hysteresis effect. 
Partial replacement of iron by Mn leads to TiFe^^Mn^ type of 
intermettalic . The hydrogen to metal ratio for x=0.7 slightly 

. * rs <*> 

exceeds 1.0. 

(c) AB z Compounds : 

These hydrides posses high hydrogen storage capacity and show 
resistance to impurities. However high thermal stability has 
resulted in their limited use. 

AB^ compounds exhibits either cubic C15 or hexagonal C14 
structure. Both of these structures have such interstices that are 
all formed by tetrahedra. As a consequence the hydrogen absorption 
increases the size of the unit cell without changing the 
structure. Main compounds of this group are ZrB^ , ZrFe^ etc. 

(d) AB a Compounds : 

Many of the intermettalic of the types AB 5 , AB z , AB fl are 
closely related to each other structurally. For example AB fl 
structure. Uhile one third of the A atoms have environment 
identical with those of A atoms of the AB 5 compounds, two third of 
the A atoms have coordination similar to that found in AB g 
systems. Therefore AB structure can be represented as 

a 9 
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AB = 1/3 (AB} +2/3 (AB ) 

3 5 2 

The AB g compounds form quit stable hydrides of the type 
AB g H^ ? and exhibit little or no hysteresis. The maximum number of 
hydrogen atoms by AB formulae united is reported in Er Co 
system. In this case a composition approaching to the theoretical 
maximum ( ErCo g H = ) is achieved. Study shows that these materials 
absorb high amount of hydrogen. 

As already pointed out earlier La-riched based intermettalic 
are very expensive because of high cost of both Ni and La. Hence 
to reduce the cost of these intermettalic attempts have been made 
to replace La and Ni with other metals. Thus Mischmetals has been 
substituted for La. Composition of Mischmetal is given in table 
1 . 3 <4> . It is relatively cheap material and does not reduce the 
advantages of LaNi badly. 

Sr 

1.5) PRESENT WORK : 

Enthalpy of absorption and desorption of metal hydrides is 
very much important, as this controls the stability of hydrides. 
In present work therefore attempts have been made to 

(1) calculate the enthalpy of formation of few binary, ternary 

and quaternary hydrides based on FeTi . Fortran computer 

programming was used to calculate the enthalpy. 

(2) study the possibility of heat pumping using metal 


hydrides . 



TABLE 1.3 


COMPOSITION OF MISCH METAL 


ELEMENT 

ATOMIC UEIGHT 

ATOMIC ( * ) 

Ut .% 




COMPOSITION 

Cerium, Ce 

140.12 

32.9 

41.9 

Lathanum, La 

138.91 

22.0 

27 . 1 

Neodymium, Nd 

144.24 

10.2 

13.0 

Praseodynium, Pr 

140.91 

- 

- 

Iron, Fe 

55.85 

34 . 6 

17.0 

Impurities 



1.0 
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CHAPTER 2 

CALCULATION OF ENTHALPY OF FORMATION OF METAL HYDRIDES 
Heat of formation of alloys is of great importance. It can be 
described in terms of a single atomic model. This relatively 
simple atomic model gives a complete account of the sign of the 
heat of formation, AH of a wide variety of alloy system. 

The basic assumption in this model is that the atomic cells 
of the two constituent metals in the alloys are to a first 
approximation identical to the atomic cells of the phase metal. 
The heat of formation alloying is then due to the change in the 
boundary conditions at the contact surface of the dissimilar cells 
in the alloy. The energy effects are described by the two terms 

(1) a negative energy term arising due to difference in chemical 

4k 

potential (A® ) between the two types of the atoms in an alloy 

(2) positive energy resulting due to the discontinuity in the 
density of electrons(An ) at the boundary between dissimilar 

V* 

1 

atomic cells. 

i 

i 

The enthalpy of formation, AH of an alloy of two transition 

I 

metals can be written as , 

I 

jj- = f (C ) g [ -ep(A® ) + Q o (An^^ ) ] 

a 

where 

N = Avogadro number 
o 

e = elementary charge 

A®* = difference in chemical potential for electron at two 


types of atom. 
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= discontinuity in the electrons density at the 
between dissimilar atomic cells. 

P = 0.147 V"‘ cm"* (du)" 1/a 
Q o = 9.4 e V* / (dul 2/a 
where 

du stands for density units (number of electrons 
volume) 

jfc 4fe 

A$ = $ ~ $ 


. i /a 1/9 

An = n 

Vfi A 

V & 


n 


1/3 


= C* C* [ 1 + 8 (C* C If ] 

A B A B 


2 ( C V 2 ' 3 4 C V 2 ^) 

A A B B 2 , 1/3 

£ s cm du 

— 1/3 -1/3 

n + n 

A B 


2/3 

c v 

A A 


_ .JZy'S _ ..2/3 

C V + C V 


A A 


B B 


c a 

B 

Where 

ion of metal A 
ion of metal B 
on of metal A 
on of metal B 


C fl = surface concentrat 

A 

C° = surface concentrat 

B 

C = atomic concentrati 

A 

C = atomic concentrati 

B 


C V 


,2 /S 


c v 2 

A A 


/a 


+ C 


v2/a 

B 


boundary 


per unit 
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V = molar volume of metal A 

A 

V = molar volume of metal B 

B 

= electron density of metal A 
n = electron density of metal B 

B 

on substituting the values of above constraints in proper 
units the expression for enthalpy reduces to 

AH =22.995 f(C°) g[-p(A4* ) 2 + Q (An 4/3 f ] KCals/mole. 

O vfi 

2.1) BINARY HYDRIDES: 

According to Miedema the metallic hydride is formed in to 
distinct steps . first, molecular hydrogen is transformed into 
metallic hydrogen, then the metallic hydrogen is alloyed with 
transition metal. 

2.2) TERNARY HYDRIDES: 

Van Mai et al extended the atomic model to the formation of 
ternary hydrides of the transition metal (AB^H^ ) on the basis of 
the following assumptions 

(a) The hydrogen atom will surround the metal atom of the 
element .which form a stable binary hydride. 

(b) The contact surface area of A-H, and B^-Hlis equal. Then 
AH of the ternary hydrides can be written as 

AHCAB H )=AH(AH )+AH(B H )-A(AB ) (ii) 

\n 2m m r. m n 

where 

AH( AH )=Enthalpy of formation of binary hydride of A 

m 

AH( B H )=Enthalpy of formation of binary hydride of B 

Vi IY* 

AH( AB )=Enthalpy of formation of original intermetallic 




117471 


compound . 
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In present work enthalpy of formation was calculated using 
Fortran programming of the following 

(i) Binary int ermetal 1 i c compounds of the type Ti^^M^where 
M=Fe/Mn/Cr/Ni/V/H and X varies from 0 to 1 with interval 

of 0.05. (Appendix a.l) 

(ii) Binary hydrides of the type M H^ where 

M=Fe/Mn/Cr/Ni/V/H and X varies from 0 to 1 with interval 
of 0.05. (Appendix a. 2) 

(iii) Ternary hydrides TiFeH and TiFeH^. (Appendix a. 3) 

(iv) Quaternary hydrides of the type TiFe^^M^H and T i F H 2 

where M=Fe /Mn/Cr /Ni /V/H and X varies from 0 to 1 with 
interval of 0 . 0 5 . (Appendix a. 4) 

Van Mal <i> found experimentally that by replacing 20% of Fe by Mn 
in TiFe, the activation procedure to form hydride could be made 
considerably simpler with reference to activation pressure etc. 
Therefore Mn can be treated as very important addition in Ti-Fe-H 
system. He also found the replacing part of Fe , in TiFe by V, Cr, 
Mn, results in stability of the hydrides. 

In present work attempt has been made to depict the effects 
of the additive metals on a purely theoretical basis. 

2.3) COMPUTATION METHOD 

The entire calculation of enthalpy of formation is 
based on the basic equation 

_ 4 n. i yg •> 

AH =22.995 f(C ) g[-p(AS ) + Q 0 On wa ) 3 KCals/mole. 


(ii) 
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* 

standard values of Q , n , and V were used in the 

vs m 

computation. These values for different metals are given in table 

< 9 > 

2.1 .Here hydrogen is considered as an imaginary metal and the 

* 

values Q =5.05 V and n =4.6 du has been assigned to hydrogen 

Vfi 

metal, which acts as the reference. 


a) BINARY INTERMETALLIC COMPOUNDS: 

Enthalpy of f ormation(AH) of binary compounds (Ti H^) was 
calculated for different M( Fe , Mn , Cr ,Ni , V , H) and X was taken from 0 
to 1 with interval of 0.05. The basic variable in this case is 
atomic concentration of the metal M. The change in AH values due 
to variation of atomic concentration is given in table 2.2 and the 
same were plotted in fig 2.1. 

b) BINARY HYDRIDES: 

AH for binary hydrides M H was calculated for M=Fe, Mn, 

l-X X 

Cr, Ni, V, Ti and X=0 to 1 with interval of 0.05. Calculated 
values are given in table 2.3 and the same are plotted in fig 2.2. 

c) TERNARY HYDRIDES: 

Enthalpy of formation of ternary hydrides TiFeH and TiFeH^UAS 
calculated using equation ( i i ). Calculated values are 

AH C TiFeH) = -10.6193 KCals/mole 
AHCTiFeH^) = -20.14448 KCals/mole 

d) QUATERNARY HYDRIDES: 

According to Van Mal a> the enthalpy of formation of quaternary 
hydrides of the typeCTiFe M H) and (TiFe, MH) can be 

i~X X X X 


expressed as 
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TABLE 2.1 


PARAMETERS FOR THE CALCULATION OF ENTHALPY OF 
FORMATION OF ALLOYS OF TRANSITION METALS 


ELEMENT 

Q* 

LS 3 

n 

vs 

(du) 

V 

m 

(cm S ) 

V 2X9 
m 

(cm 2 ) 

Ti 

3 . 65 

1 . 47 

10.5 

4.8 

Fe 

4 . 93 

1.77 

7.1 

3.7 

Mn 

4 .45 

1 .61 

7 . 4 

3.8 

Cr 

4 . 65 

1.73 

7.1 

3.7 

Ni 

5.20 

1 . 75 

6.5 

3.5 

V 

4 . 25 

1 .64 

8.3 

4.1 

H 

5.00 

1 .66 

14.3 

5.9 


22 9 

N.B For Li the value of n is l(du}, about 6 * 10 el/cm 

vs 
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TABLE 2.2 

CALCULATED ENTHALPY OF FORMATION (In KCals /mole ) FOR (Ti M 3 

i-x x 

BINARY INTERMETALLIC COMPOUNDS 

M = Fe/Mn/Cr/Ni/V/H 



X 


. 0000 
. 0500 
.1000 
.1500 
.2000 
.2500 
.3000 
.3500 
.4000 
.4500 
.5000 
.5500 
.6000 
.6500 
.7000 
.7500 
.8000 
.8500 
.9000 
.9500 
1.0000 



Ti V 

i-x x 


■15.0949 
•15 . 5069 



-12 . 4006 
-10.7918 
-9.0352 
-7 .2148 


-1.8200 

.0000 


.0000 
-.7724 
-1.5261 
■2 . 2854 
-3.0568 
•3.8306 
-4.5838 
■5.2834 
•5.8899 


•6 .7504 
•6 .6140 
-6 . 2448 
-5.6561 
•4 . 8806 
■3.9686 
■2.9815 
•1 . 9784 


.0000 

.0000 

.0000 

.0000 

-.4368 

-.3517 

-1.5245 

-.0921 

-.8631 

-.6948 

-3.0107 

- .1821 

-1.2932 

-1.0405 

-4.5049 

-.2731 

-1 . 7300 

-1.3918 

-6.0225 

- .3655 

-2.1672 

-1.7441 

-7.5503 

- .4573 

-2.5910 

-2.0870 

-9.0498 

- .5452 

-2.9820 

-2.4055 

-10.4614 

- . 6246 

-3.3175 

-2.6816 

-11.7106 

-.6906 

-3 .5741 

-2.8964 

-12.7144 

-.7385 

-3.7304 

-3.0319 

-13.3900 

-.7644 

-3.7695 

-3.0734 

-13.6639 

-.7656 

-3.6812 

-3.0113 

-13.4820 

-.7410 

-3.4645 

-2.8432 

-12.8192 

-.6913 

-3.1285 

-2.5752 

-11.6878 

- .6193 

1 — II M 


-10.1433 

-.5296 

-2 . 1855 

-1.8069 

-8.2834 

-.4279 

-1.6405 

-1.3575 

-6 . 2371 

-.3205 

-1.0886 

-.9007 

-4.1386 

- . 2128 

-.5470 


-2.0754 

-.1071 

.0000 

.0000 

.0000 

.0000 
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TABLE 2.3 

CALCULATED ENTHALPY OF FORMATION (In KCals/moIe ) FOR (H H ) 

i“X x J 

HYDRIDES 

M = Ti/Fe/Mn/Cr/Ni/V 



. 0000 
. 0500 
. 1000 
.1500 
. 2000 
.2500 
.3000 
.3500 
.4000 
.4500 
.5000 
.5500 
.6000 
.6500 
.7000 
.7500 
.8000 
.8500 
.9000 
.9500 
1.0000 


£ - x x 

i. — X X 

. 0000 

. 0000 

- 2 .2198 

.1783 

- 4.4156 

.3568 

- 6 . 6531 

.5370 

- 8.8708 

. 7076 

- 10.9480 

.8554 

- 12.7509 

. 9699 

- 14 .1615 

1 . 0451 

- 15.0949 

1.0795 

- 15 . 5069 

1 . 0752 



- 12.4006 

- 10.7918 

- 9.0352 

- 7.2148 

- 5.3922 

- 3.5976 

- 1.8200 

.0000 


1 . 0370 
.9710 
.8841 
. 7830 
. 6735 
. 5604 
.4469 
.3348 
.2242 
.1136 
. 0000 


. 0000 
-.4363 
-.8724 
- 1 .3135 
- 1 . 7339 
- 2 . 1016 
- 2.3905 
- 2.5848 
- 2.6790 
- 2 . 6770 
- 2.5892 
- 2.4305 
- 2.2175 


- 1 .6932 
- 1 . 4093 
- 1.1239 
-.8417 
-.5634 
-.2855 
.0000 


.0000 

-.1239 

-.2478 

-.3730 

-.4915 

-.5941 

-.6736 

-.7259 

-.7498 

-.7468 

-.7202 

-.6744 

-.6141 

-.5438 

-.4678 

-.3892 

-.3104 

-.2325 

-.1557 

-.0789 

.0000 


. 0000 
. 0589 
.1179 
.1772 
.2326 
.2795 
.3146 
.3365 
. 3451 
. 3414 
.3272 
. 3049 
.2765 
. 2442 
.2097 
.1744 
. 1391 
.1043 
.0699 
.0354 
.0000 


V H 

l-X X 


. 0000 
-.8828 


• 2.6552 
• 3.5194 
- 4.2949 
■ 4 . 9269 
- 5.3768 
- 5.6254 
- 5.6722 
- 5.5320 
- 5.2306 
- 4.8004 
• 4.2765 
• 3.6924 
■ 3.0776 
- 2.4544 
- 1.8366 
- 1.2282 
- .6221 
.0000 
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AHCTiFe M H) = AH(TiH ) + A H ( F e 11 H ) - A H ( T i F e K ) 

i-K X 0.5 1-x X O. 5 1-X X 


=AH(TiH ) + [ AHCFe H ) + AH(N H ) 

0 . 5 i-x O. 25 x O. 25 

- AHCFe H )] - [ AHCTi Fe ) + AH(Ti K )] 

1-K x O. 5 i-x o. 5 x 


similarly 


AHCTiF.^H H.) =AH (TiH) * [ AHCF.^H. „ ) - ABCH H, B 3 

- AHCFe^MJ] - t AH(Ti o = F Vy ) * AH(Ti 0 ,H )] 

For different M=Mn/Cr/Ni /V and K=0 to 1 with interval of 0.05 

enthalpy of quaternary hydrides were calculated. The variation in 

AH values with x are given in table 2.4, 2.5 for TiFe^^N^H and 

TiFe M H respectively and same are plotted in fig 2.3 and 2.4. 
i-x x 2 


2.4) RESULTS AND DISCUSSION 

j 

Fig 2.1 shows that the all the curves are In negative AH 
region. Hence for all the systems, formation of binary 
intermetallic compounds is possible. Stability of the various 
systems can be described in decreasing order as Ti H > Ti Ni > 
Ti-Fe > Ti-Mn > Ti-Cr > Ti-V. The stability of Ti-Ni is 

slightly 

greater than that of Ti-H system beyond the atomic fraction of 
0.6 . Atomic fraction of different for which AH is minimum 

(stability is maximum) is different. 

Fig 2.2 shows that formation of binary hydrides of Ni and Fe 
is not stable as AH is positive in most of the composition range. 
Stability of other binary hydrides in the decreasing order can be 





TABLE 2.4 


CALCULATED ENTHALPY OF FORMATION (In KCals/roole ) 


(TiF V* M * H > COMPOUNDS 


M = Mn/Cr/Ni/V 




. 0000 
.0500 
. 1000 
. 1500 
. 2000 
. 2500 
.3000 
.3500 
.4000 
.4500 
.5000 
.5500 
.6000 
.6500 
.7000 
, 7500 
,8000 
,8500 
,9000 
9500 


TiFe Mn H 
i-x x 


10 .6193 
10.6088 
10.6215 
10 . 6399 
10.6612 
10.6939 
10 . 7504 
10.8428 
10.9818 


-11 . 4336 
-11.7608 
-12.1634 


-13.2087 

-13.8509 


-15.3318 

-16.1280 

-16.9250 



-10.6193 
-10.4123 
-10.2277 
-10.0561 
-9 . 9052 
-9 . 7887 
-9 .7208 
-9 .7145 
-9 . 7809 
-9 .9290 
-10.1662 
-10.4979 
-10 .9275 
-11.4560 
-12.0812 


-13.5888 

-14.4380 

-15.3152 

-16.1918 


Ni H 

i-x x 


-10.6193 
-9 . 1380 
- 7 . 6866 
-6 . 2326 
-4.8137 
-3.4764 
-2.2588 
-1.1880 
-.2810 
.4522 
1.0072 
1.3832 
1.5829 
1.6118 
1.4788 
1 .1973 
.7865 
.2732 
-.3083 
-.9243 


TiFe V H 
*-x x 



-11.5809 

-11.8544 

-12.0915 


-12.5788 
-12.8967 
-13.3020 
-13.8139 
-14.4437 
-15 .1952 
-16.0648 
-17.0429 
-18.1135 
-19.2550 
-20.4398 


-22.8179 


1.0000 


uxcua 


•17.0791 


-1.5822 


•24.0063 
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Fig 2.3 PLOT OF AH (TiFe 
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TABLE 2.5 


CALCULATED ENTHALPY OF FORMATION (In KCals/mole ) 


(TiFe^JB^) COMPOUNDS 


M = Mn/Cr/Ni/V 



TiFe 


.0000 
.0500 
.1000 
.1500 
. 2000 
.2500 
.3000 
.3500 
.4000 
.4500 
.5000 
.5500 
, 6000 
,6500 
,7000 
,7500 
8000 
8500 
9000 
9500 


Mn H 

-X X 2 


.1448 
.1624 
.1992 
.2587 
. 3527 
.4912 
.6812 
.9273 


,6008 

,0339 

,5347 

1054 

7470 

4583 


- 26 .0635 
- 26 .9282 
- 27.8037 
- 28.6681 


TiFe Cr H 

*■ -X X 2 


- 20.1448 
- 19 . 9648 
- 19.8077 


- 19.5631 
- 19.5017 
- 19 .4987 
- 19.5652 
- 19.7100 
- 19 .9405 
- 20 .2621 
- 20.6784 
- 21.1910 
- 21.7987 
- 22.4966 
- 23.2755 
- 24.1197 
- 25.0068 
- 25.9090 
- 26.8038 



Ni H 

l-x x 2 


- 20.1448 
- 18.6898 
.2679 
- 15.8449 
- 14 . 4567 
- 13.1493 


- 10.9207 

- 10.0448 

- 9.3441 

- 8.8230 

- 8.4816 

- 8.3165 

- 8.3205 

- 8.4826 

- 8.7868 


- 9 . 7272 
- 10.3011 
- 10.9053 


TiFe 


V H 

-X X 2 


- 20.1448 

- 20.4797 

- 20.8359 

- 21.2189 

- 21.6295 

- 22.0746 

- 22.5669 

- 23.1224 


- 24.4865 
- 25.3188 
- 26.2591 
- 27.3060 
- 28.4524 
- 29 . 6854 
- 30.9864 
- 32.3314 
- 33.6917 
- 35.0374 
- 36.3494 


1.0000 


- 29.5384 


■ 27.7101 


- 11.5548 


6579 









AHCTiFe 
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described as Ti-H > V-H > Hn-H > Cr-H 

Fig 2.3 and 2.4 show that TiFe^ H H and TiFe M H systems 
partial replacement Fe by other transition metals ( Hn,Cr,Ni,V 
etc.) gives pronounced effect on the stability . 

In TiFe^^H^H system, AH is confined in negative zone for H = 
Nn , Cr , V where as for H=Ni , the enthalpy changes sign from 
negative to positive values over the composition range twice. 
Hence stability of TiFe^^H^H first decreases up to atomic 
fraction x=0 . 4 , after that it assumes positive sign up to x = 0 . 8 8 
after which it again becomes negative slowly. Stability of 
TiFe Hn H and TiFe V H first increases slowly up to x=0.-5, 
after that stability increases very fast. For TiFe Cr H stability 

A-x x 

decreases very slowly up to x=0.5 , after that it increases very 
fast. The stability of the different systems in decreasing order 
can be written as TiFe V H > TiFe Hn H > TiFe Cr H. 

l~x x A-x x A-x x 

Fig 2.4 gives the plot of TiFe H H for different x. The 

A-x x 2 

shape of the curves for different H is same as that of TiFe H H 

A-x x 

system. The only difference is that all the curves are in negative 
AH zone. Hence formation of all stable hydrides is possible. for 
H= Ni , the stability first decreases up to X=0.52 , after that it 
increases. For TiFe V H and TiFe Hn H stability increases 

A-x x St A-x x 2 

slowly up to X=0.4 after that it increases very fast. The 

stability of different systems in decreasing order can be 


described as TiFe VH > TiFe Hn H > TiFe Cr H > TiFe Ni H . 

i-x X 2 4-x X 2 4-X X 2 4-x X 2 

Hence one can see that relative stability of the different 


hydrides is same in both the system. 
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CHAPTER 3 

I 


APPLICATION OF METAL HYDRIDES IN HEAT PUMPING 
ABSTRACT: 


Hydriding of inter-metallic releases heat energy and 
dissociation of intermetallic hydrides absorbs heat energy as 
these reactions are exothermic and endothermic .The heat of 
reaction AH is related to pressure of reaction 
(absorption/desorption ) thermodynamically by the relation 
lnp=AH/RT -AS/R <4> 

Thus it is possible to utilize this reaction for various heat 
pumping application. In the present chapter the principle of heat 
pumping for cooling is discussed and calculation for 
MmNi Mn H <40> type metal hydride system developed indigenously 

5-x x <S 

done . 


The thermal energy involved in absorption and desorption and 
desorption process makes it possible with appropriate combination 
of metal -hydrogen system ,to use it at a high temperature for 
pumping heat from a low temperature source and vice-versa . This 
process can be thus used both for cooling and for heating 
applications . 

For metal hydrides the equilibrium pressure in the two phase 
region varies strongly with temperature and represented by 

lnp=AH/RT -AS/R <41 


where 

r 

AH= Heat of hydride formation (Cal/ mole of H fi ) 
p= the pressure at which hydride is formed (atm) 

AS= Chance in entropy durins hydride formation CCal /mol e°10 



36 


AS is equal for all metal -hydrogen system and constant in 

small temperature range considered and also AH is independent of 

temperature, the effect of temperature and hydriding pressure of a 

metal -hydrogen system is represented graphically by fig 3.1 <4> . 

This shows the plots of lnp^ VS 1/T for different systems. All 

2 

straight lines converge to the point Inp = -AS /Ron the line 

2 

1/T=0. This representation helps in finding the characteristics AH 
of a particular metal-hydrogen system, required for heat pumping. 

A metal hydride heat pump consists of two hydrides with 
dissimilar thermal stabilities situated in different sections of a 
sealed container, such that hydrogen gas can flow freely from one 
section to the other. Fig 3 . 2 <5> i 1 lust rat es heat pump operation in 
cooling mode. Cooling is a result of the dissociation of less 
stable hydride M H ( endothermi c reaction) at point A. The 
dissociated hydrogen flows from point A to the initially hydrogen 
depleted M at point B to form more stable hydride MH .Low 

1 ± K 

temperature hydride absorb heat from the air. By holding M^H^at 
T (i.e. maintained at ambient temperature by removing absorption 

m 

heat), P will be greater than P and hydrogen will flow from A to B 
2 £. 

and the low temperature will be maintained. Uhen the hydrogen 
content of K H has been exhausted, it is regenerated by allowing 

2 X 

it to return to temperature and contacting with H 2 after 

heating M H toT (e.g. using industrial waste heat, solar energy 

1 H H 

etc which liberates H ) at point C. Thus hydrogen will flow back 

2 

to the hydrogen depleted I^Calong C-D) since P + > P 9 -By operating 
an identical two-hydride system in tandem, continuous cooling is 
obtained. The line sketch of proposed system is given m Fig3.3 . 



AS .AH 
~R~+RT 



1/T(K~') 


T FOR DIFFERENT SYSTEMS 
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The theoretical maximum COP of the hydride based heat pump is 
given by AH /AH ‘ => 

where AH^ = enthalpy of formation of more stable hydride and 
AH a = enthalpy of formation of less stable hydride. 

In Metal Hydride heat pumping utilizing AB type of hydrides 

5 

AH^ is always lower (magnitude of negative value) than AH^ in 

order to maintain hydrogen flow from A to B (Fig 3.2 <5> ) as the 

dissociation pressure of M H is higher than that of M H . Since 

2 y ix 

the AS value for both the hydrides are about the same, this can be 
accomplished if AH^ is more negative (higher absolute value) than 
AH^ . As a result maximum value of the COP is always less than 
unity. The operating temperatures of a metal hydrides heat pump 
are dependent on the enthalpies (AH) of the hydrides chosen for 
the system. 

As already noted above, the operating temperatures of a metal 
hydride heat pump are determined by the enthalpies (AH) of the 
pair of hydrides chosen for the system. One of the great assets 
of the hydride heat pump is large range of enthalpies for 
different metal hydrides. The large number of hydrides allow the 
system to be selected to meet the heating and cooling needs of the 
user. Further, the enthalpies of ternary alloys can be varied 
over a small range by varying the composition of the alloy, 
allowing for a fine tuning of the heat pump. 

The main advantages of the hydride based heat pumps are that 
these are pollution free and can use low grade energy. By proper 
choice of two hydrides even compressor can be eliminated. 
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CHAPTER 4- 


DESIGN CALCULATIONS BASED ON Mm Ni Mn 

■*.2 o . a 


and Mm Ni 

4. 2 


INTERMETAL LICS FOR HEAT PUMPING APPLICATIONS 


A1 

o 


a 


The two indigenously developed intermetallics* 10> using Indian 
Mischmetal are MmNi Mn and MmNi A1 . Their thermodynamic 
parameters evaluated by Mangole et al <i0> are 
AH, . for MmNi Mn H = 28.40 kJ/mole of H = AH 

dofiorption -4.2 O. 3 y £ a 


AH for MmNi A1 H 

abeorption 4.2 O. 8 y 


= 32.36 kJ/mole of H = AH 


Hydrogen to metal ratio for both cases H/M =0.3 

Let us assume that absorption/desorption reactions take 10 minutes 
to complete up to 90% capacity of metal hydride. 


4.1) COST CALCULATIONS : 


Due to desorption of MmNi Mn H cooling effect per 

4 .2 O. fi y 


second, AH 


28.40 * O. P _ . ftr 

= 0.423 kU/mole of H 

<$o * iO 2 


Now for 1 ton = 3.5 kU of cooling total moles of H^ required = 


Total c oo l Lr«g 
Cool ing/mol* of H 


4 t on 

Am 7 " 


3. 5 kV 


O. 03 kv/moU 


^ 81 mole 


£ £ 

Hence metal required for 1 ton of cooling 

Bi 


Hydrogon 

0.3 


0.3 


at 270 moles 


■+ weight of hydride required for 1 ton of cooling = + ^ 


= 81 * 2 + 270 * 248.3 = 67203 gm ^ 67 kg 

Now weight of metal M required = 270 * 243.3 = 65691 gm — 66 

4 

kgAssuming coat of both the metal and metal hydride - Rs. 500/kg 
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Installation cost, pipe, etc. = Rs . 10,000 /- 

4 Total cost = Material cost + Installation cost 
= [ 500 ( 66 + 67 )] + 10000 
= Rs. 76,500 /- 
4.2) COP CALCULATIONS : 

Coefficient of Performance of the cycle is given by 
COP = AH / AH = - ? 8 - — = 0.877 

2 A 92. 3 <5 


< 5 > 
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APPENDIX 

a.l) FORTRAN PROGRAMME USED TO CALCULATE ENTHALPY OF FORMATION 

OF (Ti M ) SYSTEM 

i-x x 


print4 

4 format (78(’_')) 
x=0. 000000000001 
hh=h(x ,14.3,1.66,5.00) 
hfe=h(x,7. 1,1.77,4.93) 
hmn=h(x , 7 .4,1.61,4.45) 
hcr=h(x,7.1,1.73,4.65) 
hni=h (x , 6.5,1.75,5.2) 
hv=h(x, 8.3,1.64,4.25) 

print 6 , x , hh , hf e , hmn , her , hni , hv 

6 format ( ’ | ’ , 7 ( F8 . 4 , 2X , ’ | ’ ) ) 
x=0 . 0 

do 10 i=l , 20 
x=x+0 . 05 

hh=h(x, 14.3,1.66,5.00) 

hf e=h(x, 7.1,1.77,4.93) 

hmn=h(x, 7 .4,1.61,4.45) 

hcr=h(x,7.1,1.73,4.65) 

hni=h(x, 6. 5, 1.75,5.2) 

hv=h(x ,8.3,1.64,4.25) 

print 5 , x , hh , hf e , hmn , her ,hni , hv 

5 format ( ’ | ’ ,7(F8.4,2X, ’ | ’ )) 

10 continue 

print7 

7 format (78C’-’)) 
stop 

end 

function h(cb , vb , nb , qb) 

real h , f , cas , ebs , £ , ca , cb , va, vb , na , nb , p , qa , qb , q 

ca=l-cb 
va=10 . 5 
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na=l . 47 
qa=3 . 65 

£=2*(ca*va**.6667 + cb*vb**.6667)/(na**(-l) +nb* * (- 1 ) ) 

cal= ca* (va** . 6667 ) 

ca2= ca*va* * . 6667+cb*vb** . 6667 

cas= cal/ca2 

cbs=cb*vb** .6667*cas/ ( ca* va* *.6667 ) 
f=cas*cbs*(l+8*(cas*cbs)**2) 
p=0 .147 
q=p* 9 . 4 

h=22.9995*f*g*(-p*(qb-qa)**2+q* (nb-na) * * 2 ) 

return 

end 
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a.Z) FORTRAN PROGRAMME USED TO CALCULATE ENTHALPY OF FORMATION 

OF CM H ) SYSTEM 

!-X X 


print 8 

8 format (78( ’_') ) 
x=0 . OOOOOOOl 
hti=h(x,10.5,1.47,3.65) 
hf e=h(x, 7.1,1.77,4.93) 
hmn=h (x, 7. 4, 1.61,4. 45) 
hcr=h(x, 7 .1,1.73,4.65) 
hn i=h(x,6.5,1.75,5.2) 
hv=h(x, 8.3,1.64,4.25) 

print 7 , x.hti ,hf e ,hmn,hcr , hni ,hv 
7 format (’| f , 7 ( F8 . 4 , 2X , ’ | ’ ) ) 

x=0 . 0 

do 10 1=1,20 
x=x+0 . 05 

ht i =h (x ,10.5,1.47,3.65) 
hf e=h (x ,7.1,1.77,4.93) 
hmn=h(x , 7.4,1.61 ,4.45) 
hcr=h(x, 7.1,1.73,4.65) 
hni =h(x ,6.5,1.75,5.2) 
hv=h(x, 8.3,1.64,4.25) 
print 5 , x , ht i , hf e , hmn , her , hni , hv 
5 format ( ' | ’ , 7 ( F8 . 4 , 2X , ’ | ’ ) ) 

10 continue 

c print* ( 60 ( ' - ’ ) ) 

print 9 

9 format ( 78( ’ - ’ ) ) 
stop 

end 

function h(cb, va,na,qa) 

real h, f ,cas,cbs,g, ca , cb , va , vb , na , nb , p , qa, qb , q 

ca=l-cb 

vb=14 . 3 

nb=l . 66 

qb=5 
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fi = 2* (ca*va**. 66 6 7+cb*vb**. 6667)/ (na** (-1 )+nb** (-1 ) ) 

cal= ca* ( va* * . 6 6 67 ) 

ca2= ca*va** . 6667+cb*vb** . 6667 

cas= cal/ca2 

cbs = cb* vb* *.6667* cas / ( ca* va* *.6667) 
f=cas*cbs*(l+8*(cas*cbs)**2) 
p= 0 . 147 
q=p* 9 . 4 

h=22 . 9995*f*£*(-p*(qb-qa)**2+q*(nb-na)**2) 

return 

end 
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a. 3 > FORTRAN PROGRAMME USED TO CALCULATE ENTHALPY OF FORMATION 

OF (TiFe M H) SYSTEM 

i-K x 


print 4 

4 format ( 6 9 ( ) ) 
x=0 . 00000001 
hhmn=hh(x, 7.4,1.61,4.45) 
hhcr =hh(x , 7.1,1.73,4.65) 
hhni =hh(x , 6 .5,1.75,5.20) 
hhv=hh(x , 8.3,1.64,4.25) 
print 5 , x , hhmn , hhcr , hhni , hhv 

5 formate |’ ,5(F8.4,3X, ’ | ’)) 

x=0 . 0 

do 10 i=l , 20 
x=x+0 . 05 

hhmn=hh(x , 7.4, 1.61,4.45) 
hhcr=hh(x , 7.1, 1.73,4.65) 
hhni =hh(x ,6.5,1.75,5.20) 
hh v =hh (x, 8. 3, 1.64, 4. 25) 
pr int 6 , x , hhmn , hhcr , hhni , hhv 

6 formate | ’ ,5(F8.4,3X, ’ | ’)) 

10 continue 

pr int7 

7 f ormat ( 6 9 ( ’ - ’ ) ) 
atop 

end 

function hh(x,v,n,q) 
real n 

hl=h(l. 0,0. 5, 10. 5, 14. 3, 1.47, 1.66, 3. 65, 5. 00) 

h2=h( (l-x), 0.25, 7. 1,14. 3, 1.77, 1.66, 4. 93, 5. 00) 

h3=h( 0.5, ( 1-x) ,10.5,7.1,1.47 ,1.77,3.65,4.93) 
h4=h(x,0.25,v, 14.3,n,1.66,q,5.00) 
h5=h((l-x),x,7.1,v,1.77,n,4.93,q) 
h6=h(0.5,x,10.5,v,1.47,n,3.65,q) 
hh=hl+h2+h4-h5-h3-h6 


return 
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end 

function h(ca,cb,va,vb,na,nb,qa,qb) 

real h,£,cas,cbs,fi,ca,cb, va , vb , na , nb , p , qa , qb , q 

g=2*(ca*va** . 6667 + cb*vb** . 6667) / (na* * ( - 1 ) +nb* * (- 1 ) ) 

cal= ca* ( va* * . 6 6 6 7 ) 

ca2 = ca*va** . 6667 + cb*vb* * . 6 6 67 

cas= cal/ca2 

cbs = cb*vb* * . 6 667*cas/ (ca*va* * . 6667 ) 
£=cas*cbs*(l+8*(cas*cbs)**2) 
p=0 .147 
q=p* 9 . 4 

h= 2 2 . 9995*£*g*(-p*(qb~qa)**2+q* (nb-na) * * 2 ) 

return 

end 
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a . 4 ) FORTRAN PROGRAMME USED TO CALCULATE ENTHALPY OF FORMATION 

OF (TiFe M H ) SYSTEM 

i-x h a 


print4 

4 £ormat(69( )) 
x=0 . 00000001 

hhmn=hh (x , 7.4,1.61,4.45) 
hhcr=hh(x, 7.1,1.73,4.65) 
hhni=hh(x, 6. 5,1 .75,5.20) 
hhv=hh (x, 8. 3, 1.64, 4. 25) 
print 5 , x , hhmn , hhcr , hhni , hhv 

5 formate I’ ,5(F8.4,3X, ’ | ’)) 

x= 0 . 0 

do 10 i=l , 20 
x=x+ 0 . 05 

hhmn=hh(x , 7 .4,1.61,4.45) 
hhcr=hh(x, 7.1 ,1.73,4.65) 
hhni=hh(x , 6.5,1.75,5.20) 
hhv=hh(x, 8. 3, 1.64,4.25) 
print 6 , x , hhmn , hhcr , hhni ,hhv 

6 formate I* ,5(F8.4,3X, ’ | ')) 

10 continue 

print7 

7 £ ormat ( 69 ( ' - ’ ) ) 
atop 

end 

function hh(x,v,n,q) 
real n 

hi =h(l. 0,1. 0,10. 5, 14. 3. 1.47, 1.66, 3. 65, 5. 00) 
h2=h( ( 1-x) ,0.5,7.1,14.3,1.77,1.66,4.93,5.00) 
h3=h(0.5,(l -x) ,10.5,7.1,1.47,1.77,3.65,4.93) 
h4=h(x,0.5,v,14.3,n,1.66,q,5.00) 

h5=h((l-x),x,7.1,v,l- 77 » n * 4 - 93 ' < l^ 

h6=h(0.5,x,10.5,v,1.47,n,3.65,q) 

hh=hl +h2+h4-h5-h3-h6 

return 
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end 

function h(ca, cb, va, vb ,na , nb , qa , qb ) 

real h,£,cas, cbs , g , ca , cb , va , vb , na ,nb , p , qa , qb , q 

g=2*(ca*va** . 6667+cb*vb** . 6667 ) / (na* * (-1 )+nb** (-1 ) ) 

cal= ca* ( va* * . 6 6 6 7 ) 

ca2= ca* va* * . 66 67 + cb*vb* * . 6667 

cas= cal/ca2 

cbs = cb* vb* * . 6667 * cas / ( ca* va* * .6667) 

£=cas*cbe*(l+8* ( cas*cbs ) * * 2 ) 
p=0 .147 
q=p* 9 . 4 

h=22.9995*f*g*(-p*(qb-qa)**2+q*(nb-na)**2) 

return 

end 
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